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ABSTRACT 

We use the results of a set of three-dimensional SPH-Treecode simulations which 
model the formation and early evolution of disk galaxies, including the generation 
of heavy elements by star formation, to investigate the effects of dust absorption in 
quasar absorption line systems. 

Using a simple prescription for the production of dust, we have compared the 
column density, zinc abundance and optical depth properties of our models to the 
known properties of Damped Lyman a systems. 

We find that a significant fraction of our model galaxy disks have a higher column 
density than any observed DLA system. We are also able to show that such parts of 
the disk tend to be optically thick, implying that any background quasar would be 
obscured through much of the disk. This would produce the selection effect against 
the denser absorption systems thought to be present in observations. 
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1 INTRODUCTION 

Damped Lyman Alpha (DLA) systems are clouds of gas with 
column densities greater than 2 x 10 20 Hlcm -2 which have 
been detected in QSO absorption line systems. 

Their precise nature is still unknown, and they have 
been i nterpreted as the p r ecurso r s of present-day ma ssive 
disks dProchaska fc Wolfd Jl998h . iFerrini et alJ <1997l) L as 
sub-galactic lumps still in th e process of forming g alaxies 
through hierarchical merging (Hachnclt et al. (1998)), and 
as low-surface brig htness galaxies jjimenezetalT lll999')). 
iPettini et al] (|l997) showed that there was a large spread 
in metallicity at any given redshift for DLAs, indicating 
that they were drawn from a population of galaxies of vary- 
ing morphology and differing stages of chemical evolution. 
iLu et ail il996T> showed that iron abundances were 1/100 
solar or less for systems at z > 3, but for z < 3 many had 
abundances ten times larger than this, implying that star 
formation had begun in earnest at about z = 3. 

An argument in favour of the disk precursor idea is the 
fact that the density of gas in the damped lyman alpha ab- 
sorbers is similar to the density of matter in the form of stars 
in present-day galaxy disks. This number is approximately 
five times higher than the density of matter in the form of 
neutral hydrogen in the disks of present-day galaxies. The 
implication is that the gas in the damped Lya systems at 
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redshifts of z ~ 3 turns into stars, leaving the gas fraction 
in the z = disks at 10-20 %. 

If the damped Lya systems are massive disks at 
high redshift, then the profiles of metal absorption 
lines may be interpreted as being evidence for rotation 
]Prochaska fc W olfe (1998)). If they are smaller clumps in 
the process of merging, the velo city profiles can also b e ex- 
plained by clump bulk motions llHaehnelt et al] dl998ll ). 

In recent years, numerical simulations have shown that 
DLAs can be produced naturally wit hin the Cold Dar k 
Matter picture of str u cture formation iKataetal 
IGardner et"aH dl997l) . IGardner et"aH d200ll) L iKatz et al] 
l^9^T "showed that in a standard cluster-normalised CDM 
model (SCDM), the number of DLA systems produced 
was close to observations, although the less-dense Lyman 
Limit (LL) systems were significantly under-produced. They 
found that the lines of sight which produced damped ab- 
sorption typically passed near the centers of relatively 
dense, massive proto -galaxies with a size of ~ 10 kpc. 
IGardner et al] <ll997l) used a high-resolution SPH simula- 
tion to calibrate the relation between the absorption cross 
section and circular rotation velocity of dark matter halos 
(the (a,v c ) relation) , and then used t his w ith the Press- 
Schecter formalism JPress fc Schecterl dl974l) ) to calculate 
the number of DLAs per unit redshift in a number of differ- 
ent cosmological models. They showed that the observations 
could be fitted best by COBE-normalised CDM (CCDM) 
or SCDM, although it is well-known that CCDM produces 
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excessively massive clusters at z = 0, and that SCDM dis- 
agrees with the results of COBE. iGardner et al.l (1200 ll) car- 
ried out numerical simulations of 5 variants of the CDM 
cosmology. They found a clear relationship between the HI 
column density of each absorber and its projected distance 
to the nearest galaxy, with DLA absorption occuring for 
distances of less than 10-15 kpc, and LL absorption out to 
distances of 30 kpc. Considering only the structures resolv- 
able in the simulation, all models underproduced both DLA 
and LL absorption abundances. In order to calculate the ab- 
sorption from lower-mass halos, they calibrated the (a,v c ) 
relation using halos in the simulations and then extended 
down to lower mas ses using a mass function provided by 
Ijenkins et al.l J2001). In doing this they show that all mod- 
els are able to replicate observations of both DLA and LL 
absorption abundances so long as absorption is produced by 
halos down to v c ~ 30 — 50km s _1 . The question of the ef- 
fect of low-mass halos awaits numerical simulations which 
are able to resolve structures down to v c ~ 30km s _1 . 

In this work we too assume that DLAs are produced 
by the precursors of modern-day disk galaxies, but we re- 
investigate the possible importance of the presence of dust 
in such galaxies. Dust can increase the optical depth of the 
galaxy disk, perhaps to the point where the background 
quasar becomes obscured. We add to our analysis a sim- 
ple model for the production of dust, and we calculate the 
optical depth through a set of galaxy disks which this pro- 
duces. 

The remainder of this paper is organized as follows. In 
section 12.11 we discuss the code and in section 12.21 we dis- 
cuss the initial conditions used. In sections 12.31 and 12.41 we 
describe how we calculate number densities and metallici- 
ties along lines of sight. In sect ion 17151 we introduce our dust 
model, and in section [3] we discuss the results. Our conclu- 
sions are presented in section [I] 



2 THE MODEL 
2.1 The code 

The code which wa s used to perform the simulations is de- 
scribed in detail in IChurches et all <|200 lT) and only a brief 
overview will be given here. 

Th e code is a standard SPH-Tr e ecode iBarnes fc Hud 
<198d) . iGineold fc Monaehanl <1977ft . iLucvl (Il977ft ) which 
has been augmented to include the effects of star forma- 
tion and chemical evolution. The star formation rate per 
unit volume is given by a Schmidt Law <SchmidtHl959h L 



dpg/dt — —kp" g 



(1) 



where p g is the gas density and k and n are constants. Ob- 
servations by Kennicutt (1998) suggest that n ~ 1.5. 

Chemical evolution is modelled by an extension to the 
well-known simple model, which takes a fixed volume of the 
galaxy and turns gas into stars within it. This model makes 
the prediction that the abundance of heavy elements within 
the volume varies with time via the following relation: 
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of gas which has just been turned into stars in time inter- 
val dt, and the yield p is the ratio of the mass of metals 
which is returned to the ISM to the mass which is lost to 
the stellar graveyard, per generation of star formation. The 
extension which we have made is to let the stars and gas de- 
couple from each other and evolve as dynamically separate 
systems, whereas in the simple model they forced to stay to- 
gether in the closed box. This is achieved by applying Eqs. 
Q and to each SPH gas particle. We then form a popula- 
tion of collisionless stellar particles where those gas particles 
with a larger star formation rate are more li kely to spawn a 
new s tellar particle, as is described fully in lChurches et aD 
(2001). 



2.2 The initial Conditions 



We fo llo w the ap proac h o f iKatd jl99lTl. iKatz fe Gunnl 
<199lD. iKatzl Jl992l) . ISteinmetz fc Mulled |l994). 
Vedel et alJ (Il994fl. ISteinmetz fc Mulled il995t) . 
Williams fc Nelson! I2OO1I) and IChurches et alJ (1200 d) 
which is to model the collapse of an isolated spherical 
perturbation. These models have a resolution which is less 
than 1 kpc, and so can be used to probe internal structure, 
and have had great success in forming realistic disk galaxies. 

We start with a spherical dark matter halo which is 
composed of 6000 collisionless particles, and which is in solid 
body rotation (see below). Embedded in this are 6000 SPH 
gas particles which will model the gaseous component. As 
the simulation proceeds, gas is converted into stars, and col- 
lisionless star particles are formed. 

The effect of tidal torques by neighbouring halos on the 
dark halo to be modelled is approximated by setting the 
system in solid body rotation. The amount of rotation is 
characterised by the dimensionless spin parameter A, which 
is given by 

1/2 G' 1 M~ 5/2 



\ = J\E\ 



(3) 



where g is the total mass of gas within the box, a fraction 
Z of which is in the form of heavy elements. dS is the mass 



where J is the total angular momentum, E the total energy, 
M the total mass of the system and G the gravitational 
constant. The theoretical prediction of Peebles (1969) is that 
A ~ 0.08, a value which has been confirmed by N-body 
experiments. 

The papers of IKatz fc Gunnl <199lfl . IKatzl <1992ft . and 
ISteinmetz fc Mulled (Il995h add small-scale noise to their 
initial conditions in an attempt to model the expected CDM 
power spectrum. However, it was shown by IKatz fc Gunnl 
( 1991) that the final properties of the galaxies were essen- 
tially independent of the amount of small-scale noise. For 
this reason, we choose to describe the small-scale noise in 
the system at the start of the simulation using simple Pois- 
son noise. Experiments have been performed to verify that 
all of the physical properties of our numerical galaxies are 
independent of the particular realisation of the noise. 

The remaining parameters which specify the system 
are k and n in the Schmid t Law. Analytical work by 
lEdmunds fc Greenhowl lll995l) has shown that in order for 
the heavy element abundance to decrease with radius, n 
must be greater t han unity, while observations described by 
iKennicufra dl998l) indicate that n ~ 1.5, which is the value 
we have used. 

If one makes the assumption that a gas disk with a mass 
of M = 5 X 1O 1O M , a radius of 10 kpc and height of 200 pc 
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run 


total mass Af(xlO n M ) 


Spin Parameter A 


1 


5 


0.06 


2 


5 


0.09 


3 


5 


0.12 


4 


2.5 


0.09 


5 


10 


0.09 



Table 1. The simulation parameters. We varied the total mass 
M and the dimensionless spin parameter A. 



turns 90% of its mass into stars over a timescale of 5 Gyr 
with n = 1.5, then the value for k thus obtained is 



k = 1.45 M~ 1/2 pc 3/2 Gyr' 1 



(4) 



which we use here. 

A much wider range o f initi al conditions wer e in- 
vestigated in IChurches et alJ (120011) . IChurchesI (Il999l) and 
I Williams fc Nels^nT(l2"o"oif) . These simulations have shown 
that the parameters which have the most significant ef- 
fect upon the final properties of the model galaxies are 
the total mass of the system and its angular momentum. 
Therefore here we present results for three different val- 
ues of each. We choose three values for the total mass: 
M = 2.5 x 10 n M Q ,M = 5 x lO n M and M = 1 X 1O 12 M . 
For one of these masses (M = 5 x 10 1 M ), we varied the 
initial angular momentum, producing the following values 
of the dimensionless spin parameter: A = 0.06, A = 0.09 and 
A = 0.12. 

The simulation parameters are shown in tabled Runs 
1,2 and 3 investigate the effect of varying A while keeping M 
fixed, and runs 2, 4 and 5 investigate the effect of varying 
M for a fixed A. 



galaxy disk so that the absorption is averaged over a linear 
region in the ISM of at least a few hundred parsecs. This 
is comparable to the resolution of our simulations, so our 
column density calculations are sufficiently resolved. 



2.4 Calculation of metallicites 

Also of interest to us is the gas metallicity Z along any 
line of sight, which is defined as the ratio of the total mass 
of metals along the column to the total mass of gas along 
the column. Given that in the simulations we know both 
the mass m,i and metallicity Zi of each SPH gas particle 
(and therefore the mass of metals associated with each gas 
particle is rriiZi), this is a simple calculation. 

We need to compare our results to observations of real 
damped Lyman alpha systems. However, instead of Z, ob- 
servers usually measure the abundance of a particular heavy 
element X with respect to hydrogen, [X/H]. 

Often [Zn/H] is used as an abundance indicator be- 
cause zinc is not believed to become significantly depleted 
onto dust grains. Therefore a measurement of [Zn/H] for 
the gas is more likely to represent the true zinc abundance 
in the object. Pagel et al. (1992) give an estimate for the 
relationship between the overall heavy element abundance 
Z and O/H, Z ~ 23(0/H). Using this, we can express the 
abundance of Zn with respect to hydrogen, under the as- 
sumption that all the heavy elements increase in step with 



one another. The result is (IChurchesI il99 
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1.77 . 



(5) 



2.3 Calculation of number densities 

The surface density is defined as the mass per unit area of 
a galaxy. If one knows the mass of each of the particles in 
the galaxy, it is trivial to convert the surface density into 
a column density, which is the number of particles per unit 
area. 

Numerically we use SPH to calculate a smoothed esti- 
mate of the surface density throughout our simulated galax- 
ies and then we convert this into a column density of hydro- 
gen, using the relation 10 21 HIcm -2 = 8M pc -2 . 

For each of our models, we choose a particular time in 
its evolution. We then orient it so that it is face-on in the 
(x, y) plane, and we lay down a square grid of equally spaced 
points which cover the galaxy. The grid is a square of side 
30 kpc, centered on the galaxy. There are 16 x 16 = 256 grid 
points altogether, each separated by a distance of 2 kpc. 

We then evaluate the column density at each sample 
point, and from this we can construct the distribution of 
column densities for that galaxy. This procedure is then re- 
peated for each galaxy at different times and for different 
inclination angles. 

In real systems, the QSO may be regarded as a point 
source and the sight lines are very small in width. However, 
the sight lines pass through at least the thickness of the 



2.5 A simple dust model 

We use an elementary model for the absorp tion properties 
of dus t in g alaxies. This mo del is presented in lMathlin et alJ 
feOOll) and lChurche'sl ( 1999), and so only a summary of the 
main result is given here. 

Our model calculates the optical depth r which is pro- 
duced by looking through a column of gas which has column 
density N gas and metallicity Z. We assume that a constant 
fraction n of the heavy elements are in the form of dust, an 
assum ption for which p lausible physical arguments can be 
made fedmundsl (l200lT) 1 . The dust is uniformly mixed with 
the gas and is composed of dust particles with radius rdust 
and density pdust- This gives 

377^1.3m H N gas 

r = — 'a • W 

'iTdustPdust 

For r d ust and p d ust we have taken r d ust = 1 x 10~ 7 m 
and p dust = 2000kg m~ 3 , and for n. lEdmunds fc Ealesl 
(1998) suggest that r\ — 0.5 is acceptable. So 

r = 4.07 x 10" 20 cm 2 Z N gas (cm -2 ) (7) 

Using equation[7|we can put an estimate on the optical 
depth produced by looking through a region of a galaxy if it 
has metallicity Z and column density N ga s- This equation 
may be expressed in terms of [Zn/H] instead of Z using 
equation 
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Figure 1. Column density against zinc abundance for observed 
damped Lyman alpha systems. References to the sources of data 
given in Mathlin et al 2001. See text for details of the dotted and 
dashed lines. 
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Figure 2. Column densities and zinc abundances of lines of sight 
through a model galaxy with mass M = 5x lO n M0 and A = 0.06. 
See text for a description of this figure. 



3 RESULTS 

The simulations p resented here are draw n from a set previ- 
ously described in lChurches et all (l200lT) and a detailed de- 
scription of their evolution is given there. However, a short 
recap will be given here. 

The systems start off as uniform density spheres of gas 
and dark matter with a radius of 100 kpc. They are rotating 
about the z-axis, and are about to collapse. Collapse oc- 
curs approximately on a free-fall timescale, which is of the 
order of 1 Gyr. Material from above and below the x — y 
plane collide, the gas particles do not interpenetrate and 
so a flattened disk is formed. Rotational support about the 
z-axis gives the disk a radius of ~ 15 kpc. Knots and arm- 
like structures are formed, which subsequently merge as the 
system settles down. By this time, because of the high gas 
densities, approximately 40% of the initial gas mass has been 
turned into stars. 

As a natural result of the collapse, the gas assumes an 
exponential surface density profile, and any subsequent star 
forma tion occurs in an exponential disk. In Churches et al. 
feOOlf) we showed that these models produce galaxies which 
have negative abundance gradients in the range (—0.02 — > 
—0.07) dex/kpc, which is in broad agreement with observa- 
tions. 

We present snapshots of the absorption properties of 
our models at times of 1, 2, 3 and 5 Gyr. To give a simple 
indication of the approximate redshift reached for each of 
these times, we have used the equations for the Einstein- 
de-Sitter solution with SI = 1 and Ho = 75kms _1 Mpc _1 . 
Table [5] shows how, in this case, the time elapsed in the 
simulation is related to the redshift reached. Observations 
typically show the majority of DLAs lying in the range 1 < 
z < 4, so the overlap of the models with observations is 
primarily for the first 2 Gyr of each model. 



For each snapshot we lay down a grid which covers the 
gaseous extent of the galaxy and then we pierce the galaxy 
with 256 sight lines. For each sight line we calculate the 
column density and the zinc abundance, and from these two 
quantities we then calculate the optical depth using equation 

HI 

Figure shows the column density-abundance plot for 
the observed DLA systems. Running along the top of fig- 
ure is the gas metallicity Z, where solar composition is 
Z ~ 0.02. The upper horizontal line shows the position of 
a column density of 4.5 x 10 21 HIcm -2 , the column den- 
sity of the densest damped Lyman alpha system observed. 
The lower horizontal line corresponds to a column density 
of 2 x 10 20 HIcm -2 , which is the lower limit for classifica- 
tion as a damped Lyman alpha system. The vertical dotted 
line shows the highest zinc abundance so far observed in a 
damped Lyman alpha system, that of (Zn/H) /(Zn/ H)q ~ 
0.54. The sloping dotted curves show column density against 
abundance for constant optical depth values as expected 
from the simple dust cloud model. From the bottom up- 
wards they correspond to optical depths of 0.01, 0.1, 0.5 
and 1.0 respectively. 

One immediately notices from this diagram that no sys- 
tems lie above the r = 0.5 line. Therefore, interpreted in 
terms of our dust model, the observations all have optical 
depths less than r = 0.5. 

Figures E3 Efl ED an d ED show the results from the simu- 
lations. The top left, top right, bottom left and bottom right 
panels of these figures show the simulations at 1, 2, 3 and 
5 Gyr respectively. 

Figures|21Eil a ndEB s li ow the results of varying A for fixed 
M, and figures ED ED arl d ED snow the results of varying M for 
fixed A. 

In all cases the models fall within the box defined by 
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Figure 3. Column densities and zinc abundances of lines of sight 
through a model galaxy with mass M = 5xlO n M and A = 0.09. 
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Figure 5. Column densities and zinc abundances of lines of sight 
through a model galaxy with mass M = 2.5 x 1O U M and A = 
0.09. 
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Figure 4. Column densities and zinc abundances of lines of sight 
through a model galaxy with mass M = 5xlO n M and A = 0.12. 
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Figure 6. Column densities and zinc abundances of lines of sight 
through a model galaxy with mass M = 1 X 10 12 Mq and A = 0.09. 



the observations. This is consistent with the hypothesis that 
DLAs are indeed galaxy disks in different stages of evolution. 
At lGyr of evolution, heavy element abundances are small, 
and the sight lines occupy the left side of the diagram. As 
time increases, stars are formed, heavy elements produced, 
and the points move to the right. 

Although the models overlap with the observations, 
they occupy a greater region in the diagram. As the sight 



lines move to the right of the diagram, for a given column 
density the optical depth increases. Therefore at early times 
it is possible to have a high column density sight line which 
is optically thin due to the fact that the metallicity is small. 
However, as time proceeds the models begin to have optical 
depths of t = 0.5 and greater. The comparison between the 
models and observations can only be taken for t = 1 Gyr 
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t/Gyr M = 2.5 x 10 n M Q M = 5 x 10 n M Q M = 1 x 10 12 M C . 



4.1 
1.9 
1.2 
0.8 
0.3 



5.4 
2.2 
1.3 
0.8 
0.4 



7.1 
2.4 
1.4 
0.9 
0.4 



and t = 2 Gyr, since after this time the models have evolved 
past z = 1. By t = 2 Gyr, all of the models have a significant 
number of sight lines which have column densities greater 
than the observed upper limit. In addition, these lines lie 
above the r = 0.5 line, and so have significant optical depth. 
Very few sight lines have evolved to the maximum observed 
metallicity by this time, which may explain the paucity of 
such systems in observations. 

Figure |7| shows the fraction of sight lines which have 
an optical depth greater than 0.5 and 1.0, as a function of 
time, for all of the models. This figure shows that in all cases, 
the number of sight lines with optical depths greater than 
r = 0.5 and r = 1.0 grows sharply during the first 2 Gyr 
of evolution, and the number remains approximately con- 
stant thereafter. This is because most of the heavy element 
synthesis occurs in the first 2 Gyr of the simulation, i.e. the 
period in which the optical depth grows coincides with the 
redshift space which the observations inhabit, and it is valid 
to compare the two. 

Figure shows that the models produce galaxies which 
have ~20 - 60% of the sight lines through them with an 
optical depth greater than 0.5 and ~10 - 30% greater than 
1.0. Figure0shows that for a fixed mass, smaller A produces 
more lines with higher optical depths. This is to be expected 
since a system with less angular momentum collapses down 
to higher gas densities, and so produces more metals. This 
figure also shows that for a fixed A, larger M produces higher 
optical depths. This is again to be expected since the more 
massive systems have higher gas densities, and so produce 
more metals. 

Since the sight lines are uniformly distributed through- 
out the galaxy, and each one represents the same amount 
of surface area of the galaxy disk, the fraction of the total 
number of DLA sight lines which have a given optical depth 
is the same as the fraction of the surface area of the galaxy 
with that optical depth. Therefore, for the range of param- 
eters considered, we can produce galaxy disks which have 
20-60% of their surface area having r 0.5, and 10-30% 
with r > 1.0. 

We agree with the suggestion of iPei fc Falll <ll995f) that 
a galaxy disk with a significant optical depth could obscure 
the light from the background quasar to the point that it 
becomes unobservable. This would mean that the absorber 
would not make it into any observational DLA sample, and 
could explain why the high column density, high metallicity 
system s are not present in figure Q Calculations by Baker 
(2000) indicate that an optical depth of r = 0.5 could suffice 
to produce this selection effect. 



Table 2. The approximate redshift reached by the models as a 
function of simulation time elapsed. 
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Figure 7. Fraction of sight lines with optical depth greater than 
0.5 (left figures) and 1.0 (right figures). The top two graphs show 
the variation with A for fixed M = 5 X 10 11 Mq. Here the dotted 
line has A = 0.06, the short-dashed line has A = 0.09 and the 
long-dashed line has A = 0.12. In the bottom two graphs we 
show the variation with M for fixed A. Here the dotted line has 
M = 2.5 X 10 11 M Q , the short-dashed line has ,M = 5 X lO 11 !^ 
the A = 0.09 and the long-dashed line has M = 1 X 1O 12 M . 
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Figure 8. Fraction of sight lines with optical depth greater than 
0.5 (left figure) and 1.0 (right figure). The graphs show the varia- 
tion with inclination angle i for fixed M = 5 X lO n M0. Here the 
dotted lines have i = 0° , the short-dashed lines have i = 60° and 
the long— dashed lines have i = 80° . 



3.0.1 Inclination Effects 

To investigate the effect of inclination upon our results, we 
have repeated the above analysis for one of our models with 
two additional inclination angles of 60° and 80° . The results 
of this experiment are shown in figure |H| 

Figure|H]shows that the angle of inclination has a signif- 
icant effect upon the sight lines which have an optical depth 
t ^ 1.0, increasing their fraction from ~ 25% for i — 0° 
to ~ 35% for i — 80°. Since the distribution function of 
inclination angles is a monotonically increasing function of 
the angle of inclination, one is more likely to observe highly 
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inclined systems, and hence more likely to encounter higher 
optical depths. Therefore the figures for the face-on models 
represent lower bounds on the fraction of the disk which has 
r > 0.5 or r > 1. 



4 CONCLUSIONS 

We conclude that if the models presented here are represen- 
tative of the observed DLA systems, a significant fraction 
of the sight lines through such objects will pass through re- 
gions with column densities in excess of the observed upper 
limit of 4.5 x 10 21 Hlcm -2 . 

If we assume that half of the mass which is in the form 
of heavy metals is in the form of dust, then we predict that 
in order for the simulations to produce such a sig ht line, the 
optica l depth along that line would be r ^ 0.5. IPei fc Fall 
( 1995) have suggested that this may cause the light from 
the background quasar to be sufficiently dimmed so that 
the quasar would drop out of the observational sample, thus 
making the galaxy itself invisible in absorption. 

By putting many sight lines through different parts of 
model galaxies we have been able to calculate the fraction 
of the surface area of the disks which have an optical depth 
greater than r = 0.5 and r = 1. For our face-on models 
we find that ~ 20 - 60% has r ^ 0.5 and 10 - 30% has 
t ^ 1. If we increase the angle of inclination from 0° to 80° 
then the fraction of the disk which has r ^ 1 increases to 
~ 40%. Since highly-inclined systems are more likely to be 
observed, our face-on figures are lower bounds. 

In addition, the most optically-thick region will typi- 
cally be the inner part of the galaxy, and so the observer 
may be in danger of drawing conclusions about the proper- 
ties of these systems, based only on observations of the outer 
parts. 

In conclusion, we support the hypothesis that the ob- 
served upper limit to the column density of damped Lyman 
alpha systems is a selection effect caused by obscuration by 
dust. This implies that there could be significant numbers 
of high-metallicity, high redshift absorption systems as yet 
unobserved due to obscuration of their background QSO. 
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